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ABSTRACT 

The  spectral  properties  of  Jacobi  and  periodic  Jacobi  matrices 
art'  analyzed  and  algorithms  for  the  construction  of  Jacobi  and  periodic 
Jacobi  matrices  with  prescribed  spectra  are  presented.  Numerical  evidence 
demonstrates  that  these  algorithms  are  of  practical  utility.  These  algo- 
rithms have  been  used  in  studies  of  the  periodic  Toda  lattice,  and  might 
also  be  used  in  studies  of  inverse  eigenvalue  problems  for  Sturm-Louivi lie 
equations  and  Hill’s  equation. 
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Signi f icance  and  Explanation 


In  this  report  we  present  algoritlims  which  solve  two  inverse  eigenvalue 
problems  that  arise  in  matrix  theory.  Computational  evidence  is  presented 
that  demonstrates  that  these  algorithms  are  of  practical  utility. 

The  first  inverse  eigenvalue  problem  considers  what  additional  information 
uniquely  determines  the  entries  of  a Jacobi  matrix  if  we  know  its  eigenvalues. 
Recall  that  a Jacobi  matrix  is  a real,  symmetric  tridiagonal  matrix  whose  next 
to  diagonal  entries  are  positive. 

The  second  inverse  eigenvalue  problem  considers  what  additional  information 
uniquely  determines  the  entries  of  a periodic  Jacobi  matrix  if  we  know  its 
eigenvalues.  A periodic  Jacobi  matrix  is  obtained  by  replacing  the  entries  in 
the  upper  right  and  lower  left  corners  of  a Jacobi  matrix  by  the  same  positive 
number . 

Inverse  eigenvalue  problems  of  this  nature  arise  in  mathematical  physics. 
For  example,  the  construction  of  a linear  array  of  masses  interconnected  by 
springs  with  presci rbed  normal  modes  of  vibration  leads  to  such  inverse  eigen- 
value problems.  In  addition,  the  construction  of  a ladder  network  of  inductors 
md  capacitors  with  prescribed  transmission  characteristics  also  leads  to  such 
inverse  eigenvalue  problems. 

Finally,  FORTRAN  subroutines  which  implement  these  algorithms  are  presented 


in  an  appendix. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summaty 
lies  with  MRO,  and  not  with  the  author  of  this  report . 
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THE  CONSTKUCTION  OF  JACOBI  AND  PERIODIC  JACOBI  MATRICES 
WITH  PRESCRIBED  SPECTRA 

Warren  E.  Ferguson,  Jr. 


I.  Introduction.  A periodic  Jacobi  matrix  is  any  real,  symmetric  matrix 
of  the  form 


I,  = 


N-l 


where  b.  ->0  Vi. 
1 


This  paper  shows  how  one  can  construct  a periodic  Jacobi  matrix  with  prescribed 
spectra.  For  example,  there  is  a family  of  periodic  Jacobi  matrices  with  Vj, 


• • • , \N  as  eigenvalues  if  and  only  if  the  numbers 
be  ordered  so  that 


XL--.,*  are  real  and  can 


\. 


- s 


4—5 


Similar  problems  have  been  studied  by  other  authors  [2,  12). 

The  results  presented  in  this  paper  are  basts!  upon  an  analysis  of  the  spectral 
properties  of  periodic  Jacobi  matrices.  The  main  tool  in  this  analysis  is  the 
knowledge  of  the  spectral  properties  of  Jacobi  matrices.  Recall  that  a Jacobi 
matrix  is  any  real,  symmetric  tridiaqonal  matrix  whose  next  to  diagonal  entries  are 
positive.  Our  cannonical  Jacobi  matrix  will  be  the  matrix  obtained  by  deleting 
from  L the  last  row  and  column,  that  is 


N-2 


where  b.  ■»  0 V i, 

l 


An  algorithm  which  constructs  a Jacobi  matrix  with  prescribed  spectra  is 
presented  in  Theorem  2.  This  algorithm  is  derived  from  the  fact  that  any  real , 
symmetric  matrix  has  real  eigenvalues  and  a corresponding  full  set  of  real , 
orthonormal  eigenvectors.  We  hasten  t o point  out  that  essentially  the  same 
algor  it  let  was  presents!  by  de  Boor  and  Golub  (l).  Similar  problems  have  been 
studied  by  other  authors  [H,  0] . 


ffpohaored  by  Ihc  Jkrlfty  un3er  ConVriol  Tfol  D*AG29-75-c-6dl4.  mUi  tal 
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The  spectral  properties  of  periodic  Jacobi  matrices  are  considered  in 
Section  3.  The  results  presented  in  this  section  are  derived  from  a matrix 
analog  of  Floquet  theory  [11].  In  Section  4 wo  use  these  results  to  charac- 
terize the  family  of  periodic  Jacobi  matrices  with  prescribed  spectra. 

We  present  the  results  of  several  numerical  experiments  in  Section  5. 

These  results  demonstrate  that  the  alqoritlims  presented  in  Theorems  2 and  <3 
.re  of  practical  utility.  Indeed,  these  alqoritlims  have  been  used  in  performing 
numerical  experiments  on  the  periodic  Toda  lattice  [4].  In  Section  6 we  con- 
clude the  paper  with  several  ctmments. 


2.  Spectral  Properties  of  Jacobi  Matrices.  In  this  section  we  will  con- 
sider the  spectral  properties  of  the  Jacobi  matrix 


Observe  that  J is  a real,  symmetric  matrix.  Consequently  J has  real  eigen- 
values Ui,-*-,uN_i  and  a corresponding  set  Y^,*'*,YN_1  of  real,  orthonormal 
eigenvectors  [13,14],  If  Y denotes  the  matrix  whose  ith.  column  is  Yj  then 
Y ’ s an  orthogonal  matrix  ..nd 


JY  = YD  where 


Many  important  relationships  between  the  eigenvalues  and  eigenvectors  of  J 
can  be  derived  from  the  representation 


(ill  - J) 


-1  = Y(ul-D)"1 


(2) 


of  the  resolvent  of  J . For  example,  by  comparing  the  entries  in  row  1, 
column  N-l  of  (2)  we  arrive  at  the  identity 


r N-l 

’n-2  “ ' 1 


<>'(li) 

U“VI. 


Y Y 
l,k  N- 


1,K 


Here 


(3. a) 


u(li)  = dot  (uT-J)  (3.b) 

is  the  characteristic  polynomial  of  J and  Yj,  j denotes  the  entry  of  Y in 
row  i , column  j.  Another  important  identity,  used  by  Stieltjes  in  his  treat- 
ment of  inverse  eigenvalue  problems,  can  be  derived  from  (2)  by  comparing  the 
entries  in  row  1,  column  1 (or  row  N-l,  column  N-l.) 


t 


In  the  work  that  follows  we  will  demonstrate  that  the  entries  of  J can 
be  recovered  from  the  entries  on  the  diagonal  of  D and  in  the  first  row  of  Y . 
Before  we  describe  this  process  let  us  introduce  the  following: 


Definition  1:  (a)  The  Jacobi  matrix  J is  characterized  by  the  data 

{p,y}  if  and  only  if 

(1)  u,  , •••,p  , are  the  eigenvalues  of  J , and 

1 N-l 

(2)  y, , are  the  first  components  of  a set 

1 N-l 

Y,  , • • • ,Y.  , of  real,  orthonormal  eigenvectors  of 

1 N-l 

J corresponding  to  ^ . 

(b)  The  data  {p,y}  is  compatible  if  and  only  if 


(1) 

V 

• • • , u 

'mN-1 

are  real , 

distinct  numbers,  and 

(2) 

V 

• • • , v 
'yN-l 

are  real. 

nonzero  numbers  whose  squares 

sum 

to  one. 

We  feel  justified  in  using  the  words  "characterize"  and  "compatible"  in  this 
manner  because  the  following  theorem  is  true. 

Theorem  2:  Data  characterizing  a Jacobi  matrix  is  compatible.  Furthermore 

each  set  of  compatible  data  {p,y}  characterizes  a unique  Jacobi  matrix  J . 

The  entries  (a,b)  of  this  Jacobi  matrix  are  computed  by  the  algorithm: 


1. 

Yo,j 

= 0 V j 

2. 

Vi 

= Yj  V j 

3. 

For 

i = 1 , • • • ,N-2 

4. 

VN_1  v2 

ai  £1  UkYi,k 

5. 

b2  = ^-1[(p  -a.)Y.  - b.  Y.  ]2 

i k i i,k  l-l  i-lfk 

k 

6. 

Y.  - . = t“[  (p . -a . ) Y . . - b.  .Y.  . .] 

1+1,3  b^  3 l 1,3  l-l  i-l,3 

7. 

Next 

i 

8. 

aN-l 

vN-1  2 

^1  UkYN-l,k  * 

Proof:  The  proof  of  this  theorem  will  be  presented  as  a sequence  of  three 

lemmas.  . 


-3- 


1 twma  J.l:  bat  a characterizing  a Jacobi  matrix  is  ci<m[>at  i 1*  1 *• . 

i’lvH't  : let  tin-  Jacobi  mat  i ix  J bo  characterized  by  the  data  (u,y|.  The 

•• 1 ■ ‘tti'  nrivr.s.ii  ily  teal  Ivi’aust'  they  at  o the  eigenvalues  of  a real,  symmet t ic 
matrix,  by  definition  the  y's  at o real,  and  their  squares  sum  to  one  because 
they  may  lx-  considered  t o he  the  entries  in  t lie  first  row  of  the  orthogonal  matrix 
Y in  (11.  Consider  the  limiting  form  of  the  identity  ( t)  as  u tends  to  ti  . 

l- j were  a repeated  eigenvalue  then  ..i'(uj)  0 and  so  we  would  be  tot.  ed  t.> 

conclude  that  N****Yl-.I  0 . which  is  imi'osslble  beeause  each  b,  N 0 . rherefore 
t h«*  ii's  are  distinct  and,  as  vi  tends  to  , we  infer  that 


1 


V,  - Y.Vi.i 


Consequently  the  y's  ate  nonzero  because 


Y , 


V 1 

Y 


(41 


l.j  * 


lemma  J . J : v.iven  cimiwttihle  data  (li,yl  the  algot  itlim  ot  Theorem  ovxnput 

the  entries  (a.b)  of  a Jacobi  matrix  J ohai act er ized  by  the  data  (u,yl. 


1'rcH'f : First,  we  infer  that  this  algor  it  ten  computes  the  entries  (a,b)  ot 

s.me  Jacobi  mat  i ix  J only  it  the  value  of  b computed  in  step  S is  novel  zeto. 


from  the  oompat  ibi  l it  y ot  the  data  we  infet  that  b 
' fot  some  f ' N-l  then  step  t'  implies 


txit  b» 


Ol 


Y . 

i , 1 I 


i . 1 


<-l 


- u 


Li.  J 


Y 


I 


but  this  is  im]H>ssible,  fot  no  mat  t ix  of  ot  der 
e igenva lues . 


It  bj 


V 1 


•N-l  • d 


' . 1 J 

N-l  has  N-l  d i s t i net 


Second,  we  will  dtmonst  i at  e that  the  numbers 
r i t lim  t orm  the  entties  ot  an  ottluvtonal  matrix  Y 
satisfy  the  ot  thonormal  ity  relations 


VN-1 
} . Y 


t , k ' i , k 


l . 1 


t ot 


Y;  i com  put  ed  by  this  algo- 
, that  is  the  tow;,  ot  Y 


1 


1 , 


• , t 


(SI 


and  i - l , • • • , 

N-l  . 

From  the  c 

ixnpat  i bi  l i t y ot 

t he  dat  a 

( u , y ) 

we  inter 

t hat 

(‘<1  is  true  fot 

i 

1.  If 

(M 

is  1 1 ui>  for  i 

!,•••,< 

t hen 

t he  tv'll  ow  i 1st 

argument  d<mon; 

;t  rat 

es  t hat 

it  i 

s also  true  for 

i ♦ 1 . 

Cleat 

ly  steps 

■'  and  t 

imply  that  (‘>1 

is  t 

rue  fot 

i “ 

» * l . For  i v 

t’  step  t< 

imp! i e 

s that 

yN-1 

it1 

' <»  l 

. k'  i , k 

1 

b 

l V"-1  H V.  Y. 

j^l  k t , k t , 

v “ A * 'S  • 

i k v v , \ 

- b 

v - 

r\’  - 1 , i 1 

-4- 


» 


Tin-  i iaht  side  of  this  equality  is  zero  1 ot  1 • i because  step  4 was  execute.), 
and  it  1-.  zero  toi  t - C because  step  t>  implies  t liat 


V 


N-l 


1 kY  ( , kY  i ,k  ^1  Yi',kU'i-lY)-l  ,k  * *'  | Y j , k * V’l+l.k'  * V , M 1 


d,  wc  will  demonstrate  that  the  data  iu.yl  characterizes  .1  . It  wi  1 


be  sufficient  to  pi-'ve  that  the  matrices  J , Y constructed  by  this  alijot  i t tim 
satisfy  (1).  Step  t>  implies  that  JY  “ YP  if  we  can  show  that  the  numbers 


V,  (l  i - Vi’Vi.j  - V.'V.vi  v 1 


.it.'  zero.  Tin-  techniques  presented  in  the  previous  paragraph  can  be  used  t c ,..,n 
on  t rat e t hat 


t-N-l 


V ' ' y Y 0 

; l N.k  j.k 


t or  i 1 , • • • ,N - 1 . 


Since  the  rows  of  Y form  a real,  orthonormal  basis,  we  infer  that 


Y 0 

N.  1 


V i 


I .ewma  . ' : bach  set  of  compatible  data  characterizes  at  most  oik  ,’aoohi 

mat i ix . 


1'ioot  : let  d be  any  Jacobi  matrix  characterized  by  tin-  com pat  ible  date. 
ii.,yl.  Then  y f***,y  ^ are  the  first  components  of  a set  Y^,***,Y^  ^ oi 
real,  orthonormal  eigenvectors  of  J zori  espondinu  to  t tie  eigenvalues  j , • • • , 
j . If  Y denotes  the  mat  i ix  whose  itli  column  is  Y,  t hen  Y is  an  " t ho 


■tonal  mat  i ix  and 


J Y Y V 


where 


I 0 


N-l 


We  will  now  prove  that  the  entries  (a,b)  of  J art'  ident  ical  to  the  outlie. 
(a,b)  v't  t ht>  Jact'bi  matt  ix  J computed  by  the  uKiot  itlim  presented  in  Thc'iem 


Tilt'  ent  i ies  Y 


i , 1 


't  Y satisfy  the  orthonormal  ity  t elation; 


i-N-1 


' i , k'  i , k 


i . 1 


v i,i 


because  YY  1 . The  entries  Y.  of  Y also  satisfv  tht'  t ecut  t en  'el. 

l , 1 


b . Y . . * a . Y , * b Y u Y . 

t - 1 t - 1 , t it,)  t i ♦ 1 , 1 t i , t 


V l , i 


whet  e 


Y 0 ami 

0,  l 


Y„  • 0 

N,  1 


v i 


because  J Y Y 1'.  When  the  recur t enee  relation  is  multiplied  by  Yj 
result  is  -.limned  ovet  i we  tint!,  usino  tin'  t't  t lu'notmal  it  v t elutions,  t 


a VN_l  ; Y* 
i kl  k t,k 


L 


V I 


The  i «vm  i i‘n»’c  rol.it  ion  also  implies  that 

Y . - - [ ( l;  . -a  ) Y , - b Y . . ) 

1*1.  J b 1 1 l.j  l-l  » - 1 . J 

atul,  when  this  identity  is  squared  and  t ho  result  is  s minut'd  ovei  l , tho  ortho- 
normality  rolat  ions;  imply  that 

V^’ 1 Ml.  - a . ) Y - b.  Y 1*  . 

i : l k i i , K i- 1 i-l , k 

Starting  with  tho  fact  that 


Y,  y V i 

l . 1 l 

it  i . easily  shown  by  iiuluct  ion,  following  tho  sequence  of  c.  input  at  ions  piesented 
in  t ho  algor  it  tin,  that  tho  ontrios  (a, Id  of  j are  identical  to  tho  ontiios 
(a , b)  of  j . 

■ 

t.  Specti.il  ii  opoi  t IQS  of  Periodic  Jacob i Mat  i ices.  In  this  sect  ion  wo 
will  coniidts  tho  spivtial  properties  of  the  periodic  Jacobi  matrix 


wht'i  o h s 0 
1 


V 1 . 


Throughout  this  soot  ion  wo  will  uso  J t o i opi  esont  the  Jaeobi  mat  i i\  obtaiinxi 
by  delet  ing  1 1 om  I tho  last  row  ami  column. 

observe  that  1 is  a real,  symmetric  matrix.  Consequent ly  1 has  real 
eigenvalues  and  correspond ing  set  ot  real,  orthonormal  eigenveotois  (It,  141. 
let  r be  an  eigenvector  of  I corresjw'nding  to  t he  eigenvalue  \ . Then  the 
opponent z of  ; form  a nontrivial  solut  ion  ot  t he  recurrence  lolation 

d',  ‘ 

0 N 


b .*  » a . t . t b . / . 

i-li-l  ii  11*1 


Vz  . 


V i 


which,  satisfies  the  boundary  conditions 
■ „ and 


N 


Nt  l 


1 


by  anaK'gy  with  Kloquet  theory,  which  analyzes  the  analaious  problem  foi  ordinal \ 
differential  equat  ions  [111,  let  us  cons  idol  the  nontrivial  solut  ions  ot  the  room 
fence  lelat ion  which  sat istv  the  boundary  oondit  ions 


and 


1 


l 


Here  the  parameter  o is  called  the  Kloquet  multiplier  ot  r . This  problem  ha- 
only  the  trivial  solvit  ion  when  e 0 , while  toi  c > o a nontrivial  solut  ion 
exists  if  and  only  it  \ is  an  eigenvalue  ot  the  matrix 


With  tho-o  tact-  in  mind  lot  us  mt  uxlua'  t ho  tollownui: 

' o t 011 1 i on  <:  lot  bo  ohar.ict  or  i zed  by  t ho  vl.it  a \i,y‘  and  I avo  ..(i.l 

a i t characteristic  polynomial,  rhen  the  Floquet  multipliers  P1»****Pn  \ of 

! i’oi  t o i-v'tni  m.j  to  ■ 1 ■ * * * > i\_  i ai  e t ho  numbers  defined  by  the  J elation 


b,  •••to  u>*  (y  ' h y 

1 N i l N't 


V i . 


1’heoi  or\  -1 


■ha«  act  ei  1st  ie  polynomial  ot  l admits  the  i opt  o-ont  at  ion 


•f  l \ 1 -1  ' 1 


- • • • b (,\(\)  - (.«+  -)  ) 

IN 


wheie  \(\),  .-.tiled  the  discriminant  of  1 , is  independent  of  o . The  Floquet 

multipliei  . t , • • • , .’>j  _ j of  l,  cot  t espon.l  inn  to  the  eigenvalues  lij  # • • • #Vijj_  j 
t at  i --t  \ t he  i el  at  ion 


A ( ' . * 

l 1 


v i . 


' .«>  .•  i nenva  . nos 


\ o.t  1 ate  teal  and  can  N-  ot  doted  so  that 


S ■ S • S'  S'  s •••  • 

I’toot  -inn  olomentaty  propel t i os  ot  determinants  it  is  not  hatd  to  danoi 
t i it  o t hat 

def(\l-l  ) -b  • • • b^,  (1  - — ) . 

shot;  both  ide  at«'  mt  t\|t  at  t'd  with  respect  t o we  t i rid  that 

dot  ( \ I -1  1 h • • • b • \(  \'  - (,>♦';  1 1 . 

't  com  so  the  constant  i>t  inttxirat  ion  b •••  b^  \(1'  is  neeessai  ily  indept'ndor. 


lot  bo  ihat.u'ti't  u’txl  by  the  data  ; t , y ' . Then  i « ‘ * * < V\’- v are  the 

tu  t component of  a sot  Yj  , • • • .Yjj-i  ot  loal  , otthononn.il  oi.ionvoct  oi  s ot 
. 't  t«  ond  itui  to  it-  o iai'liva  1 lies  1. | , • • • , njj_  \ . lot  Yj  j denote  the  ith  com 
. orient  ot  V,  . ! t ont  the  definition  (t>>  of  the  Floquet  multipliois  and  t he  iden- 
tity i4'  wo  i’lt  ot  t h.at 


c. 


i'ons«\juont  lv  . , if.  an  oiuonvaluo  ot  I fot  imcIi  j and  wo  inf«‘i  from  ( ' 
t hat  (8)  i :•  t ruo  . 

■ ; i nn  t hi'  ii  t init  iv'ti  ( tO  ot  I hr  K l oquot  mu  1 1 i j > 1 i cf  s w«  li tniuc  1'  t hat 

1 (u  1.  0 V l . 

i I 

Whon  t ho  I'Mi'HV.ilui’  ■ ot  J ait-  oiiii’i  id  so  t hat 


"l  “N-l 

t In’ll  wu  uiti'i  ^81  that 


(-1)  A ( Vi  ) 

1 — 


bocauso  t ho  maututudo  of  . + - is  novor  loss  than  two.  Consoquont ly  t ho  < ‘ioon- 
valuos  \ j , • • • , of  I.,  whioft  aro  t ho  root:'-  of  N(0  aro  roal  and  fan  t'i 

ordorod  so  t hat 

'l  ' S * \>  " ••• 


- V 


boouuso  tho  ooi't  t ioiont  (l>  1 of  in  .\(\)  is  positivo.  . 

I N 

A typioal  disoi  im  inant  tot  a pot  iodio  docobi  matt  is  t ot  ordi-t  N n t 
i 1 lust  rat  i'd  in  Kiquro  1.  In  this  t iuuro  wo  dopiot  tho  rolat  ionship  hotwot-n  tho 
oiqonvaluos  \ j , • • • , \jj  i'f  I.  and  tho  Kloquot  multipliers  ^ » • • • > ojj-'j  of  1 . 

rosf'ond  inu  to  the  eigenvalues  '•  1 > ’ * * • > n-  1 of  i1  • 

hot  us  vnt  roduoo  tho  follow!  tut: 

Dof  in  it  ion  : (at  1'ho  pot  iodio  Jaoobi  matt  is  I is  oharaotoi  i;:od  hy  t he 
data  l if  and  only  it 


(1) 

A 

a 

i * **• 

* < an  • 

(.') 

B 

i •**  l> 

• t 

•i 

( n 

v 

• • 

’ * UN - l 

aro  tho  eigenvalues  ol  , and 

(41 

*V 

• • 

’ ’ ‘N-l 

aro  ttio  Klvxjuot  multipliers  of 

V 

• • 

’ ,UN-1 

Tho 

dat  a 

{ A , B , i , . 

')  is.  ov'mpat  thlo  it  and  only  if 

in 

A 

\ :» 

a roal 

numl'ot  , 

B 

i s 

a roal 

, posit  ivo  ntimbot  , 

( M 

lll  ' 

• • 

’ ’ ' N-l 

at o toal,  dist  inot  numhors,  and 

(4) 

■'l  ’ 

• • 

’,0N-l 

aro  ti'al  numbers  whioh  satisfy 

w 1 1 

h 

(i  !•••().-•,  1 . 

1 N-l 
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M<  feel  iu>.t  if  ied  in  vising  the  words  "characterize"  ami  "compatible"  ir>  '■  -i 
Runnt'i  1«  dusr  Mu-  f Mowing  theortm  is  true. 

in  oi  in  c:  I 'at  a eharact  er  iz  mq  a periodic  Jacobi  matrix  is  compat  ibl  e . 

Km  t ht  more , e.i  -ii  set  of  compatible  data  tA,B,i.,,)  characterizes  a unique  pet’, 
die  Jacobi  mat i vx  ! . The  entries  (u,b)  of  this  periodic  Jacobi  matrix  are 

computed  by  the  algorithm: 


1 . b* 


yN-1  P 

K1  V'V 


b / . ) 

N 1 1 


Recover  J from  the  data 


4.  b. 


...  h \ 

1 N-2  N 


5*  A - (ai  + •••  + Vi1 

with  (;  ) 1 ••  • v . ) . 

I N-  I 

Proof:  rhe  proof  of  this  theorem  will  be  presented  as  a sequence 

three  lemmas. 

■ 

l.emma  e.l:  I'ata  char  act  er  i z ing  a periodic  Jacobi  matrix  is  compa'  . 

Vtoot:  l.et  the  periodic  Jacobi  matrix  I,  be  characterized  by  the  iita 

A,H,.  ,.  . Clearly  A is  a real  number  because  it  is  a sum  of  real  numbs  , 

while  B is  a real,  positive  number  because  it  is  a product  of  real,  pod: iv. 
numbers.  The  . 's  ati  real,  distinct  numbers  because  they  are  the  eiqenvi'.  in 
ct  the  Jaoobi  matrix  J . while  the  definition  ( f>  ) of  the  c ' s makes 
obvious  t hat  they  are  real,  nonzero  numbers  which  satisfy  . 1 . - 0 f " 

because  ( ) is  also  the  characteristic  polynomial  of  J . 


Lemma  t . : liven  compatible  data  {A,P,n,o}  the  alqoritlim  of  Tbi  ororu  e 

•ompute  the  entries  (a , o)  i'f  a periodic  Jacobi  matrix  I characterize  b\  : in 

data  iA,P,n,ph 

Proof : The  data  li.,yl  used  in  step  .1  is  compatible,  therefore  it  i . . 1 . : 

that  this  alqoritlim  computes  the  entries  (u,b)  i'f  some  periodic  Jacobi  tin  • ■ \ 
let  I be  charact  er  i zed  by  the  data  {A,B,u,p}  • Krom  steps  -1  and  5 it  \ -’laj 
that  A A and  B B.  In  view  of  Theorem  .1  we  know  that  J is  chara  • o.  -\ 

the  data  .V'  . Therefore  n ^ = ui  for  all  i anil  from  the  definition  m t h- 
Lloquet  mult ipliers  we  know  that 


B -j  ' (i  ) l\  y:  V i . 

) l N 3 


Step  therefore  implies  that 


o . f or  all  i . 
l 1 


Lemma  i> . 1 : Kach  set  of  compatible  data  characterizes  at  most,  one  periodic 

Jacobi  matrix. 

Proof:  bet  L be  any  periodic  Jacobi  matrix  characterized  by  the  compatible 

data  {A,B,u»p)  • Let  the  Jacobi  matrix  J , obtained  by  deleting  from  L the 
last  row  and  column,  be  characterized  by  the  data  (ii,y)  . Without  loss  of  gener- 
ality we  may  assume  that  each  y-  is  positive,  for  if  y-  is  the  first  component 

J J 

of  an  eigenvector  Y:  of  J then  -y-  is  the  first  component  of  the  eigenvector 
-Yj  of  J . We  will  now  prove  that  tne  entries  (a,b)  of  L are  identical  to 
the  entries  (a,b)  of  the  periodic  Jacobi  matrix  L constructed  by  the  algor it’im 
of  Theorem  6 . 


By  definition  the  I- loquet  multipliers  P^'”*',0N-1  L corresponding  to 

u i , • • • , t<N_i  satisfy  the  relationship 

-2  -2 

B = -Pjio’  (n  . ) bN  y V j . 


The  sum  of  the  squares  of  the  y's  equals  one  because  the  data  { m , y 1 is  compat- 
ible, therefore 


~2  vN-1  B 

b = - i , 7- 

N kl  pka.  (uk> 


and 


J-  /-  2 

bN  / 


V j • 


In  view  of  steps  1 and  2 we  infer  that  bjj  = bjj  and  yj  = yj  for  all  j.  Since 
both  J and  J are  characterized  by  the  same  data  then  Theorem  2 impl ies  that 
J = J.  Finally,  steps  4 and  5 imply  that  b^-i  = and  = aN  . 


4.  Periodic  Jacobi  Matrices  with  Prescribed  Spectra.  With  these  basic- 
facts  established  let  us  now  consider  how  we  can  characterize  the  family  of  perio- 
dic Jacobi  matrices  whose  eigenvalues  are  X^,***,X  . 

Let  I,  be  a periodic  Jacobi  matrix  characterized  by  the  data  {A,B,;i,p}. 

Then  X ^ , • • • , XN  are  the  eigenvalues  of  L if  and  only  if  the  discriminant  A(\) 
of  L admits  the  representation 

A ( X ) = 2 4 ( X - X ) * * * (X  — X ) . 

B 1 N 

Therefore  the  problem  of  characterizing  the  family  of  periodic  Jacobi  matrices 
with  proscribed  spectra  is  intimately  related  to  the  problem  of  characterizing  the 
family  of  periodic  Jacobi  matrices  with  prescribed  discriminant.  Let  us  introduce 
the  following: 

Definition  7:  For  each  polynomial  p(X)  let  F(p)  denote  the  family  of 

periodic  Jacobi  matrices  whose  discriminant  is  p(X)  . 

The  problem  of  characterizing  which  periodic  Jacobi  matrices  belonq  to  F(p)  is 
answered  in  the  following: 
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N . The  dat  ,i 


t A , H , , , . 


c 


Theorem  8:  Let  p(\)  be  a polynomial  of  deqn 

characterizes  a member  of  F(p)  if  and  only  if: 


(1) 

tin1  data  {A,B,u 

, p 1 is 

compat i ble , 

(2) 

P(\)  ^ l\N  - A 

xN_1  + 

lower  powers  of  V 1 , 

and 

(3) 

p(w.)  ■ p.  ■*  — 

. Pj 

V j 

Kurt hermore , 

F(p)  is  nonempty 

i f and 

only  if 

(4) 

the  coefficient 

„ N 

of  X 

in  p(  X)  is  posit ive, 

and 

(5) 

p(,\)  has  local 

ext  rema 

at  N-  1 real  , dist  me 

t point 

\)  N . • « \* 

1 N-l 

with 

(-1)  ’ p ( v ) 1 2 V j 

Proof:  The  proof  of  this  theorem  will  be  presented  as  a sequence  of  two 

lemmas . 


Lemma  B.l:  Tin-  data  {A,B,u,p}  characterizes  a number  of  f (p)  it  and 

only  if  conditions  (1,2,3)  of  Theorem  8 are  satisfied. 

Proof:  If  the  data  (A,B,u,pl  characterizes  a member  of  F(p)  then 

Theorems  4 and  f>  demonstrate  that  conditions  (1,2,3)  of  Theorem  B are  satisfied. 

Let  us  now  suppose  that  conditions  (1,2,3)  of  Theorem  8 are  sat isf ied.  Let 
,\(A)  be  the  discriminant  of  the  periodic  Jacobi  matrix  charactei ized  by  the  data 
I A , 8 , [i , p I . Now 


q ( A ) A ( X)  - p(X) 

N-l  N 

is  a polynomial  of  degree  N-2  because  the  coefficients  of  \ , V X ( v > , 

p(\)  agree.  Theorem  4 also  implies  that  q(p.)  0 V j and  so 

q(X)  0 

because  I In'  only  polynomial  of  deqree  N-J  which  is  zero  at  N-l  dist  iiu't  p 'int- 
is the  trivial  polynomial.  Consequent  ly  the  data  {A,b,ii,pl  eharactci  i.c  a 
member  of  F(p)  • 


Lemma  B.2:  F (p)  is  nonempty  if  and  only  if  conditions  (4  ,5)  of  Theorem  H 

are  s.itisfied. 

Proof:  If  f (p)  is  nonempty  then  Lemma  B.l  and  the  mean-value  theorem  r,i 

h<-  used  to  demonstrate  that  conditions;  (4,r>)  ot  Theorem  H are  sat  i .tied 


Let  us  now  suppose  that  condit  ions  (4,5)  ot  Theot  em  8 at  e sat  ist  ied.  Let  A,b 
be  determined  so  th.it 

p(\)  — ( \N  - A \N  * t lower  powers  ot  \ I 
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c 


arul  pL 


P..  , be  solutions  of 
N-l 


Then 
dat  a 


p ( v . ) 
J 


V j • 


the  data  {A,B,v,p}  is  compatible  and 
{A,B,v,pl  characterizes  a member  of 


from  Lemma  8.1  we  infer  that 

F (p)  . 


the 


■ 


Using  Theorem  8 it  is  not  hard  to  prove  the  following: 


Corollary  The  periodic  Jacobi  matrix  I,  has  A ,•••,>  as  its  eigen- 

values  if  and  only  if 


u r (a  ) 

B '0  B 


where  Aj.(X)  2 + — ( • (X-X  ) . Furthermore,  there  is  a periodic  Jacobi 
matrix  with  X , • • • , \N  as  its  eigenvalues  if  and  only  if  the  numbers  Xj,»**,Xn 
can  be  ordered  so  that 


X2  - S 


X4  - X5 


5.  Numerical  Experiments.  Let  us  now  present  the  results  of  several 
numerical  experiments.  These  experiments  were  carried  out  on  a UNIVAC  1110  in 
single  precision  floating  point  arithmetic  (27  bit  mantissa)  using  FORTRAN  ver- 
sions of  the  algorithms  presented  in  Theorems  2 and  t>  . 

In  the  first  experiment  we  test  the  algorithm  presented  in  Theorem  2.  The 
results  of  this  experiment  are  presented  in  Table  1.  Observe  that  this  a Igor  it  tan 
has  difficulty  in  recovering  the  Jacobi  matrix  described  in  Kxamplc  '. 

Kxper iment  1: 


1.  Select  a Jacobi  matrix  J of  order  N-l. 

Compute  the  data  !p,y}  characterizing  J [11,  14,  15]: 

(a)  use  bisection  to  compute  the  p’s  , and 

(b)  use  inverse  iteration  to  compute  the  y's  . 

1.  Use  the  algoritlim  presented  in  Theorem  2 to  reconstruct  the  Jacobi 
matrix  J character ized  by  the  data  lp,y}. 

4.  Output  the  error  II J - j II  where 


In  the  second  experiment  we  test  the  algoritlim  presented  in  Theorem  6.  The 
results  of  this  experiment  are  presented  in  Table  2.  Observe  that  the  Jacobi 
matrices  used  in  the  examples  of  Experiment  1 are  obtained  by  deleting  the  last 
row  and  column  from  the  periodic  Jacobi  matrix  used  in  the  corrosj'onding  examplt  : 
of  Experiment  2. 


N . 


Fxpej  iitirnt 


1.  Select  a periodic  Jacobi  matrix  L of  order 

2.  Compute  the  data  (A,B,u,pl  characterizing  L : 

(a)  use  the  obvious  sum  to  compute  A , 

(b)  i:;e  the  obvious  product  to  compute  B , 

(c)  compute  the  data  { jj  , y > characterizing  J as  described  in 
Step  2 of  Experiment  1,  and 

(d)  compute  the  p's  using  Equation  (f>)  . 

1.  Use  the  algorithm  presented  in  Theorem  6 to  reconstruct  the  periodic 
Jacobi  matrix  I.  characterized  by  the  data  lA,B,\i,p}  . 

4.  Output  the  error  Ilh-Lll  where 

II All  = max  la  . . I . 
i,i  13 

In  both  of  these  experiments  we  have  not  worked  with  matrices  of  order 
N to.  The  reason  why  we  have  not  worked  with  matrices  of  order  N > to  may 
l«c  explained  as  follows.  In  Example  2 of  Experiment  2 some  of  the  components,  oi 
y m the  data  lu,yl  become  smaller  as  N increases.  For  example,  the  small- 
est component  of  y changes  from  2xlo_y  for  N = 15  to  2x10”*°  for  N to. 
Since  the  Eloqurt  multipliers  p depend  on  the  squares  of  the  data  y we  w'll 
run  inti'  underflow  problems  when  N •*  30  . The  immediate  remedy  for  this  under- 
flow problem  is  the  use  of  logarithms  in  the  computation  of  the  Floquct  mult  •- 
pliers.  However , underflow  also  occurs  in  the  computation  of  y when  N 5'>, 
consequently  the  use  of  logarithms  is  not  a panacea. 


t>.  Comment s . Let  m(n)  be  the  characteristic  polynomial  of  the  Jacobi 
matrix  J obtained  from  J by  deleting  the  first  row  and  column.  By  come.u  inn 
the  ■ nfr  ies  ot  (.’)  in  row  1,  column  1 we  find  that 


,e(|i) 


n 


N - 1 ii'(p ) 


II-  U. 


l,k 


k " "k 

Thi  s ident  it  y was  used  by  Stielt  jos  in  his  study  of 
As  ii  tends  to  |i  we  deduce  that 


dll.) 

I 


(o' (ll  J Y*  , 

.1  i r J 


V j 


inverse  eigenvalue  pioblom:  . 


From  this-  identity  we  infer  that  the  eigenvalues  of  J strictly  interl.ii.  those 
of  J . Furthermore,  from  the  eigenvalues  of  J and  J we  can  recovei  In  d.C  i 
lii,y)  character izing  J and  hence  J itself. 


It  is  interesting  to  note  that  the  algorithm  presented  in  Theorem  is  » 'J 
in  some  versions  of  the  implicit  shift  QR  algorithm  [ 1 t 1 . These  version:  o' 
the  ytt  alqoritlim  make  use  of  flu1  fact  that  if  B = y A y" , where  B i:  .m  unto 
duced  upper  llessenberg  matrix  and  Q is  a unitary  matrix,  then  the  entri  • o'  I 

and  i.'  are  uniquely  determined  from  the  entries  of  A and  the  entries  in  tin- 
first  row  ot  £ . In  our  application  A = D,  B = J and  Q * Y. 
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We  can  also  recover  the  Jacobi  matrix  J from  the  eigenvalues  and  the  last 
components  of  the  corresponding  real,  orthonormal  eigenvectors  of  J . To  under- 
stand why  let  us  consider  the  permutation  matrix 


0 1 


We  find  that  S*  = I , therefore  from  Equation  (1)  we  deduce  that 
(SJS) (SY)  = (SY) D . 

Consequently  the  algorithm  presented  in  Theorem  2 states  that  the  entries  of 


can  be  recovered  from  the  entries  of  D and  the  entries  in  the  first  row  of  SY, 
that  is  the  last  row  of  Y . 


It  also  appears  that  Theorem  2 can  be  extended  to  some  class  of  band 
matrices.  For  example,  let  the  real,  symmetric  matrix 


have  as  eigenvalues  and  Y^,*'*,Y>j_ i as  the  corresponding  set  of 

real,  orthonormal  eigenvectors.  If  Y denotes  the  matrix  whose  jiii  column  is 
Yj  then  Y is  an  orthogonal  matrix  and 


K Y = Y D 


where 


D 


0 


t'Ui'wiin  t hi  .ligament  j.,  , ,.nted  m 1 iira.i  .’.A  we  . » » t i » .> » .m  .ilgott*  mi  v.  1 . J 

i • i'v.i  f tim  tin-  •■tit  i if  in  P .m.i  in  t hr  f 1 1 ■ . t two  tow:  of  Y . 

P 1 ■ • .i|  • i Py  . lu!'  .nnl  Wi'l  oh  1 7 1 outlines  how  out  .in  modi  t tin  usual 

a Uh'i  1 t Ian  and  omput  ■ 1 1 t ..  t 1 \ t lie  il.it  a ; , , y I ehn t ast  et  i ;•  i nil  a . i. loot  i mat  i i x . > . 

I'iii  t lii  im.  t r,  t he  i ! pd|  et  also  present  a mat  t i x Vet  ••  i on  ol  t he  eel  el  i a t e.l  . .el  I ami  - 

! entail  ■.. . l tit  i on  to  the  t nvet  se  eigenvalue  | t oPl  ,in  I ol  a e l a • . i St  u t m-l.  i . mv  i 1 1 ■ 

i 1 '1  l,jn  . 

1'1'e  ; a[  et  l'-,  K an'i't.  let  lli'l  desctilies  an  a lgoi  i t lmi  th.it  ean  pi  in  e.l  to  ion. tin. 
a 1 1 .et  i itu  mint  wtiose  "shape"  is  [reset  i P(«d . By  tin  "shape"  ot  a .1  i se  t iin  i mint  w. 
ate  tetett  i tig  to  the  value  of  the  d i s.  i im inant  at  oaeh  ot  its  N-l  teal,  (list  in  t 
l oi  i 1 ext  tana.  pot  a p|  l i eat  i oils  ot  K.  unmet  ei  1 a 1 got  1 1 lim  to  t he  pet  toil  ie  To.  la 
lattje.  we  t et  .a  the  i ea.let  to  the  forthcoming  papet  I -1 1 . 

> f u l ml  ormat  ion  eoneet  n ing  ptopet  t ie:;  ot  pet  i . *1  t e .laeo'  j mat  l ire:,  is  eon 

t a t ne.i  in  111.  We  would  a 1 so  like  to  state  t hat  t he  ana  1 y s i pt  . sent  ml  in  Sect  t on 

< an  !e  extended  in  the  same  generality  to  "ant  i -pet  i .xi  i e"  JaooPi  mat  t tee  .!  t h< 


whet . h 


N-l  ' N 


Aokttow  1 edgemen  t ..  The  allt  hot  would  like  to  thank  I’tofessot  (’ . de  Hoot  , 
I*  • 'li  lika,  ...  doluP  and  p.  M.laughlin  lot  iwtal  i nt  ormat  ive  discussions. 
Indeed,  the  te  ults  pt  esen  t ed  in  this  papet  arose  t t om  wot  k done  with  Ptasehka 
. " l s.tgh  1 i n 1-11  on  t he  |'et  t.xlic  Tola  1 a t t i e.  while  t lie  a 1 got  1 1 Inn  | g . • on  t . -d  in 

‘•ote'u  i essentially  the  am.-  algojittam  present ixl  Py  de  Hoot  and  doluP  Ml. 
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Error 


4'l0 

2'U) 
5' 1 0 
2'10 
2 '10 
O' 10 


-8 

-7 

-7 

-7 

-7 

-7 


(N+l-l)/N-2  1 = 1 , • • • ,N-1 
1-  (N-l)  /N  1 * 1 , • • • ,N-2 


Error 

-8 

4 ' 1 0 

-7 

1 ' 1 0 
4 ' 1 0 

-7 

*'10 

-7 

*'10 

-7 

O' 10 


1 : 


Kt'sult;-.  ot  Expor imont  1 


Example  1: 


Example  2: 


A (I)  = 

-2  I = 1 , • • • 

B ( I)  = 

1 I = 1, • • • 

A(N)  = 

0 

B(N-l) 

= B(N)  =1 

N 

Error 

5 

9*10-8 

-7 

10 

5x10 

-6 

15 

1x10 

20 

2xl0~6 

25 

3x10 

“6 

30 

5x10 

A (I)  = 

(N+l-U/N-2 

B (I)  = 

1-(N-I)/N 

A (N)  = 

0 

B(N-l) 

= B(N)  =1 

N 

Error 

5 

lxio"7 

10 

2xl0~7 

15 

4xl0~7 

20 

3xlO~7 

25 

6xlo"7 

30 

4xlo~7 

,N-1 

,N-2 


I 

I 


1 , • • • ,N-1 
1 ,* • • ,N-2 


Table  2 - Results  of  Experiment  2 


r 

p 
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Example  3: 


A ( I ) 

a 

I/N-2 

I = 1, 

• • • ,N-1 

B(I) 

a 

1 - I/N  1=1, 

• • • ,N-2 

A (N ) 

s 

0 

B(N- 

l) 

- B (N) 

= 1 

N 

Error 

5 

4"lo"8 

10 

mo"7 

15 

mo"3 

20 

o 

2*10 

25 

4vl0° 

30 

5x10° 

Table  2.  Results  of  Experiment  2. 
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8.  Appendix 


In  this  section  we  present  listings  of  several  FORTRAN  subroutines  that 
the  author  used  while  performing  various  computational  experiments. 

No  warranties,  expressed  or  implied,  are  made  by  the  author  that  this 
program  is  free  of  error.  It  should  not  be  relied  on  as  the  sole  basis  to 
solve  a problem  whose  incorrect  solution  could  result  in  injury  to  person 
or  property.  If  the  program  is  employed  in  such  a manner,  it  is  at  the  user's 
own  risk  and  the  author  disclaims  al  ; liability  for  such  misuse. 
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